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Recently, new strongly interacting phases have been uncovered in mesoscopic systems with chaotic 
scattering at the boundaries by two of the present authors and R. Shankar. This analysis is reliable 
when the dimensionless conductance of the system is large, and is nonperturbative in both disorder 
and interactions. The new phases are the mesoscopic analogue of spontaneous distortions of the 
Fermi surface induced by interactions in bulk systems and can occur in any Fermi liquid channel 
with angular momentum m. Here we show that the phase with m even has a diamagnetic persistent 
current (seen experimentally but mysterious theoretically), while that with m odd can be driven 
through a transition which spontaneously breaks time-reversal symmetry by increasing the coupling 
to dissipative leads. 
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The interplay of disorder and interactions is a rich 
source of unexplained phenomena in the bulk, especially 
in two dimensions [1], despite three decades of theoreti- 
cal effort [2-4]. In mesoscopic systems one is confronted 
with phenomena not seen in the bulk, such as Coulomb 
Blockade oscillations [5,6] of the zero-bias conductance, 
or persistent currents in mesoscopic rings in a small ex- 
ternal magnetic field [7-14] (for reviews see [15]). It has 
been realized in the last several years that one needs 
to take interactions seriously in order to understand the 
experimental Coulomb Blockade peak spacing statistics 
[6]. This understanding has led to Universal Hamilto- 
nian [16,17] treatments for weak-coupling (electron gas 
parameter r s = l/aoy/np ss 1, where ao is the Bohr ra- 
dius and p is the electron density), in which a constant 
charging interaction and a constant exchange interaction 
are kept in addition to the single-particle energy. 

Recently, two of us investigated [18] the stability of 
the Universal Hamiltonian to other interactions of the 
Landau Fermi- liquid type [19], which are expected to be 
present in ballistic quantum dots with chaotic boundary 
scattering, but are not in dots deep in the diffusive limit 
[20]. The Landau interaction u(9(k) — 0(k')) depends 
only on the angle between the momenta k, k' of the in- 
teracting particles [19], and can be parameterized in two 
dimensions by a set of Landau parameters u m : 

oo 

u(9-6') = u + u m cosm(6-6') (1) 

m— 1 

The values of the Landau parameters depend on the 
strength and range of interactions, and can only be ac- 
cessed numerically [21] for large r s (strong interactions). 
It was found [18] that while the Universal Hamiltonian 
was stable in the regime u m > u* m = —1/2 In 2, it be- 
came unstable to a mesoscopic Pomeranchuk [22] phase 
for u m < u* . Soon afterwards one of us and R. Shankar 



[23] showed how to construct the large- N theory of the 
strong-coupling regime (the dimensionless conductance g 
plays the role of N). Details of this treatment will appear 
soon [24]. 

In this Letter, treating Coulomb Blockade and persis- 
tent currents within the same approach, we show that 
the mesoscopic Pomeranchuk phases display unexpected 
signatures in the persistent current, including a diamag- 
netic persistent current (seen experimentally [12-14] but 
so far unexplained) in a model without superconductiv- 
ity for m even, and spontaneous time-reversal symmetry 
breaking for m odd. 

We want the effective Hamiltonian in an energy win- 
dow of width the Thouless energy Et around the Fermi 
energy. This is the regime of validity [15] of Random 
Matrix Theory (RMT) [25]. For ballistic structures 
Et = hvp/L, where vp is the Fermi velocity and L is 
the linear system size. The dimensionless conductance is 
defined as the number of single-particle energy states (of 
mean spacing A) in this window g — Et /A. Our effective 
Hamiltonian [18,23] has a nonintcracting part represent- 
ing the chaotic scattering at the walls, and a Fermi-liquid- 
like interacting part which conserves momentum. The 
order parameter [23] cr of the Pomeranchuk phase is a 
two-dimensional vector whose magnitude a and direction 
X represent the size and direction of the maximum Fermi 
surface distortion. The shape of the deformed Fermi sur- 
face is given by acos(m9 — x), where m is the angular 
momentum channel in which the instability occurs. To 
determine the behaviour of cr the fermions are integrated 
out and an effective action is obtained [23] , the dominant 
part of which is self- averaging for large g. 

Seff = 9 2 J dt - ^) + A(<V) + 

gf^\a(u>)\ 2 f(u) (2) 
The "kinetic" term f(uS) behaves like ui 2 for ui <C A 
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FIG. 1. The effective potential in the Mexican Hat in units 
of A plotted as a function of the angle \ f° r the cases m even 
and odd for q = 20 and u m =1.7. 



FIG. 2. The ensemble-averaged minimum of the effective 
potential in units of A (for g = 20 and u m = 1.7) as a function 
of the crossover parameter <j> for m — 1 and m = 2. 



and like |cj| for w> A, indicating the Landau damping 
of cr. For large g the saddle-point of this action domi- 
nates [23], since fluctuations are down by 1/g. A strong 
enough attractive Landau parameter u m < u* m leads to 
symmetry-breaking [18,23]. 

We have numerically evaluated the effective potential 
for the collective variable cr and analyzed its dependence 
on external magnetic flux. This effective potential can be 
obtained as the ground state energy £(er) of a noninter- 
acting fermion Hamiltonian where cr = cr(cosx* +sin%j) 
appears as a parameter: 



{H a (x))ai3 = £ a S a (3 - g<7M af3 (x) 



(3) 



where the first term encodes the chaotic scattering (with 
the eigenvalues e a controlled by RMT), and M a p(x) 
represents the coupling between the collective mode cr 
and particle- hole excitations of the fermions [23]. In the 
large-g limit, the ground state energy of the system in 
the strong-coupling regime is [23] simply the value of 
the effective potential at its global minimum. This en- 
ergy automatically contains the contributions which are 
specific to the particular disorder realization in addition 
to the self-averaging contributions determnined earlier 
[23] . The averaged energy landscape of cr in the strong- 
coupling phase would be a "Mexican Hat", with rota- 
tional symmetry [23]. The sample-specific contributions 
break this symmetry completely for m even, leading to a 
single minimum, and to a two-fold symmetry for m odd, 
as can be seen in Figure 1. 

To see the relation of the above to the persistent cur- 
rent we note that [15] 



dF(<p) 
dcf> ' 



(4) 



where <fi is the external flux piercing the sample. At zero 
temperature the Free energy F is just the ground state 
energy £, so we desire to obtain £ as a function of (j). 
We obtain this by taking the noninteracting part of the 
Hamiltonian from the RMT ensemble of crossover Hamil- 
tonians [25] parametrized as 



1 + C 2 (f> 2 



H s + C^H, 



(5) 



where H$,a represent symmetric (T preserving) and an- 
tisymmetric (T breaking) random matrices drawn from 
their respective normalized ensembles [25], and C is a 
factor of order unity which depends on the shape of the 
sample and the precise nature of the chaotic scattering 
at the boundary [17]. 

Figure 2 shows the dependence of £ on the crossover 
parameter C<p for m even and odd, and clearly shows 
the diamagnetic behaviour for m even. To see why this 
is special, consider what is known about persistent cur- 
rents [15]. In a mesoscopic ring penetrated by a flux, the 
ground state energy has to be periodic in the flux, since 
an integer number of flux quanta can be gauged away. 
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dF 

((/)) = -—= I lS m(2n 



o) + h sin(47n/>/(/>o) + 



(6) 



where <f>Q — h/e is the flux quantum. Only the even 
moments l2 n survive disorder-averaging [15]. 

In the noninteracting case the typical, fluctuating val- 
ues of the Fourier coefficients are (for small n) I n .typ ~ 
Et/4>o while the average is (l2n) ~ A/</>o- Experiments 
typically measure (l2n) and a few other low-order har- 
monics [9-14]. 

Interactions, when included in renormalized first-order 
perturbation theory [8], produce (I2) ~ M*"f^" where fi* 
(of order 1) is the dimensionless Cooper-channel inter- 
action at low energies. Thus, interactions enhance the 
average persistent current, but if fi* > (I2) should be 
paramagnetic, while if fx* < it should be diamagnetic. 
This prediction [8] , while of the same order-of-magnitudc 
as the experiments [9-14], has the wrong sign. Materials 
that show no sign of superconductivity (implying that 
/j* > 0) show [11-14] a diamagnetic (12). Many expla- 
nations have been proposed to account for this puzzle 
(a recent one being Ref. [26]), but the question remains 
open, as summarized in Ref. [27]. In this context a dia- 
magnetic persistent current of order Et/4>o m a model 
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without superconductivity is striking. The fact that our 
treatment is nonperturbative in the interactions [23] en- 
ables us to evade the usual sign [8]. Our approach, while 
suggestive, is not directly applicable to the experiments 
on Au and Ag rings [12,14] since the samples are not 
likely to be in the strong-coupling regime, and are not 
fully in the ballistic limit (the elastic mean free path is of 
the same order as the system size). On the other hand, 
our theory would apply directly to ensembles of ballis- 
tic GaAs rings of the type used in Refs. [11,13], but at 
stronger coupling. 

Let us now turn to m odd. The exact degeneracy of 
the two global minima separated by 7r in the angle \ can 
be seen from Figure 1, and can be proved analytically 
using the relation H* (x) = H a (x + n) . The two degen- 
erate minima are related by the time-reversal transfor- 
mation T. A particular value of x leads to a distortion 
of the Fermi surface along the direction specified by x- 
Under T, k — * k and a distortion of the Fermi sur- 
face for odd m maps to an inequivalent state at x + 77 
with the same energy, since the underlying Hamiltonian 
is T-invariant. The ground state of a Hamiltonian quan- 
tum system with a two-fold degenerate potential is the 
symmetric combination of the two minima. This applies 
to the isolated mesoscopic structure, whose dynamics is 
Hamiltonian at energy scales smaller than A (for energy 
scales in the range A « w « Et the dynamics is dissipa- 
tive with ohmic dissipation, see Eq. (2)). The splitting 
between the symmetric and antisymmetric combinations 
is the tunneling amplitude between the two minima, here 
Ae~ 9 . The two minima correspond to states carrying op- 
posite persistent currents, and are macroscopically dis- 
tinguishable. 

The coupling of the mesoscopic structure to the leads 
produces ohmic dissipation at arbitrarily low energies. 
This is precisely the case of the Caldeira-Leggett model 
[28] considered and solved by Chakravarty [29] , and Bray 
and Moore [29] . The effective action of our model at low 
energies (ui <C A) is 



gjdt{V{ X {t)) + \{d X /dtf) 



ttA 2 



f dtdt 
J (t-t') 2 



dtdt' gin 2 / \_ : \__ 



(7) 



where X is the angle of cr in the Mexican Hat, V(x) 
is the doubly-degenerate realization-specific potential of 
Figure 1, and Y is the level width induced by coupling to 
the leads. The long-range interaction in imaginary time 
comes from a "kinetic" term, which in turn arises from 
the Landau damping of cr due to decay into particle-hole 
pairs at arbitrarily low energies, possible because each 
formerly sharp level a is broadened by coupling to the 
leads. 

The model has a weak-dissipation phase in which the 
ground state is still the symmetric superposition of the 
two minima, and a strong-dissipation phase in which the 




2.0 



FIG. 3. The ensemble-averaged second Fourier coefficient 
of the persistent current in units of A/0o as a function of 
coupling \Um\ > |«m| — l/ln2 for m — 1 and m = 2 and 
9 = 20. 



particle is localized in one minimum. The transition be- 
tween the two phases [29] occurs for g(T / A) 2 (cr 2 ) w 1. 
For large enough g even a weak coupling to the leads 
(r « A/^/g, since (cr 2 ) w 1) is sufficient to meet this cri- 
terion, and leads to localization in one minimum of the 
twofold degenerate effective potential even at zero tem- 
perature, corresponding to a spontaneous breaking of T. 
The T-breaking transition could be monitored by mea- 
suring the peak- height statistics [30] . 

If one turns on an external flux, one minimum moves 
up in energy as the flux increases while the other moves 
down. The ground state (which moves down) displays a 
paramagnetic persistent current of order Et/4>q- In the 
isolated dot, or in the case with weak dissipation, one 
starts with a symmetric superposition of the two minima 
as the zero-flux ground state. As <j) increases the system 
crosses over to fully T-broken dynamics when the energy 
difference of the two minima is greater than their splitting 
in zero field, which is Ae~ 9 . Thus, the crossover will 
occur for an external flux <j>x ~ 4>o e ~ 6 '■> as compared to 
the noninteracting crossover flux <f> w (j^o/y/g- For strong 
enough dissipation, the ground state already breaks T, 
and the variation of £ contains a term first-order in (f>, 
which implies a spontaneous persistent current at zero 
flux. 

Figure 3 shows the ensemble-averaged second Fourier 
coefficient (I2) for m even and odd. Ideally one would 
find a periodic behaviour of £ with (f>. This cannot be 
captured by the crossover Hamiltonian, but must be put 
in by hand. Since we do not know the number C con- 
necting the crossover parameter to the flux </>, there is 
an inherent ambiguity in this procedure. We show (I2) 
for three different choices for C . As can be seen, the 
qualitative results are unaffected. 

Finite temperature produces additional interesting ef- 
fects for m odd. For temperatures exceeding the tun- 
nel splitting Ae -9 the superposition states are irrelevant, 
and one can think of the two minima as separately ther- 
mally occupied. Because of their (near) degeneracy, there 
will be a Curie-like susceptibility of l/T of the persis- 
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tent current (and therefore the magnetization) to exter- 
nal flux. 

In summary, we have shown that some surprising sig- 
natures of the mesoscopic Pomeranchuk regimes show up 
in the persistent current. There is a diamagnctic per- 
sistent current (for to even) without any superconduc- 
tivity. The m odd case undergoes a spontaneous time- 
reversal symmetry-breaking transition as the coupling to 
the leads is increased, and displays a spontaneous per- 
sistent current (at zero flux) in the T-broken phase. It 
would be very interesting to explore the behaviour of 
persistent currents in the quantum critical regime [31], 
where large fluctuations of the order parameter <r and 
finite quasiparticle lifetime [32] at low energies are ex- 
pected [23]. 

We thank R. Shankar for illuminating discussions and 
E. Mucciolo for pointing out Ref. [21], and the National 
Science Foundation for partial support under grants 
DMR 98-04983 (DH and HM) and DMR 0071611 (GM). 
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